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\ Clustering is one of the most complex phenomena known to the structure of atomic 

Mh ■ nuclei. A comprehensive description of this ubiquitous phenomenon goes beyond standard 

QQ ■ shell model and cluster model frameworks. We argue that clustering is a consequence of 

■ an openness of the nuclear many-body system. To illustrate this point, we study the near- 
threshold behavior of exceptional points in ^''Ne and ^*S. 

§1. Introduction 

^ . The nuclear many-body problem is among the most difficult challenges in physics 

\ due to its paramount complexity and importance at various scales: from sub-nucleon 

to astronomical. This complexity is the reason that the theory of the atomic nucleus 
is still very much a "work in progress." A comprehensive, unified description of all 
J> \ nuclei at the nucleonic level would require merging structural approaches like ab 

■ initio, configuration- interaction (shell model), and nuclear density functional theory 
. with reaction approaches such as the continuum shell model or the modern theory 
\ of optical potential. In this quest, we need new theoretical concepts and new experi- 

fSj ■ mental data, especially close to the drip lines where non-local effects appear already 

. in the ground state. 

Since such a coherent description is not yet within our grasp, however, one av- 
enue is to create a patchwork of advanced models describing selected nuclear features. 
Consequently, our understanding of nuclear structure is incomplete and even inco- 
^ \ herent in certain aspects. Low-energy excitations in light atomic nuclei provide the 

■ case in point. Both bound states and resonances can be modeled using the nuclear 
shell model (SM) and the cluster model (CM), but the picture of how nucleons are 
organized within many-body states is different in these approaches. The degrees 
of freedom of SM are the valence protons and neutrons moving in selected shells 
whereas the CM is based on a molecular viewpoint that employs more composite 
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effective building block^ such as a— particles. Moreover, both models assume the 
nucleus to be a closed quantum system (CQS) that is completely isolated from the 
space of scattering and decay channels. Can the comprehensive understanding of 
low-energy excitations possibly emerge from such disjointed physical pictures? 

In this paper, we shall highlight some of the theoretical issues underlying the 
clustering phenomenon by analyzing properties of the continuum shell model (CSM) 
wave function^'® in the neighborhood of the proton decay threshold. Salient 
features of near-threshold CSM states will be illustrated by studying exceptional 
points-^ (EPs) of the complex-extended CSM Hamiltonian. Exceptional points are 
singularities of eigenvalues and eigenvectors for complex values of selected model 
parameters; they can be associated with square root branch point singularities of 
the eigenvalues in the vicinity of level repulsion. By investigating the energy- and 
(complex) interaction-dependence of exceptional points around the proton emission 
threshold and identifying salient features of eigenfunctions due to a coupling to a 
common decay threshold, we hope to reveal generic features of the clusterization 
mechanism. 

Our paper is organized as follows. In Sec. [2] we examine physical arguments that 
necessitate the use of an OQS description of the clustering phenomenon and reiterate 
its universal features. A brief summary of the real-energy CSM and EPs is presented 
in Sec. [3l The properties of exceptional threads, i.e., trajectories of EPs in the space 
of system energy and continuum-coupling interaction, is discussed in Sec. [H As 
an illustrative example, we study the J'^ = 0^ CSM eigenfunctions of ""^^Ne and ^^S, 
which are mixed due to the coupling to a common i = proton decay channel. These 
results provide a good insight into the configuration mixing mechanism operating in 
more complicated decay channels involving composite charged particles such as ^H, 
^He, ^He, and ®Be. Finally, Sec. [5] contains the main conclusions of this work. 

§2. Salient configuration mixing phenomena near the dissociation 

threshold 

What can be said about the structure of many-body states in the narrow range 
of energies around the reaction threshold? Are those properties very dependent on a 
particular realization of the Hamiltonian? Is there a connection between the branch 
point singularity at the particle emission threshold and the appearance of cluster 
states? In the context of these questions, the observation by Ikeda et al.'^ offers 
important insight: Q-cluster states can be found in the proximity of a-particle decay 
thresholds. Curiously, this observation has not raised much interest in studies of the 
coupling to particle-decay channels as the mechanism behind enhanced correlations 
in near-threshold SM states. 

The conjecture of Ikeda et al.-' can be formulated more generally: the coupling 
to a nearby particle/cluster decay channel induces particle/cluster correlations in 
SM wave functions which constitute the imprint of this channel. In other words, the 
clustering is the generic near-threshold phenomenon in OQS that does not originate 
from any particular property of the Hamiltonian or some symmetry of the nuclear 
many-body problem. We claim that this conjecture holds for all kinds of cluster 
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Fig. 1. Example of two Borromean systems: ^^C and ^^Li. The near-threshold states in these nuclei 
- the Hoyle resonance in ^'^C and the weakly-bound 3/2^ ground state in ^^Li - exhibit a 
strong imprint of the nearby decay threshold, i.e., ^^C— >*Be-)-a and ^^Li-^- 2n+^Li. 



states including unstable systems like dineutron, or ^Be. Figure [T] shows two exam- 
ples of the Borromean systems for which the lowest-energy threshold corresponds to 
the emission of a cluster resonance: ^Be in ^^C and dineutron in ^^Li. 

Universality of the clustering phenomenon stems from basic properties of the 
scattering matrix in a multichannel system.'^' Indeed, the decay threshold is a 
branching point of the scattering matrix. For energies below the lowest particle- 
decay threshold at energy Ei, one finds the analytic phase with a single solution 
which is regular in the entire space. Above Ei and below the energy E2 of the next 
decay channel, the new analytic phase corresponds to two regular solutions of the 
scattering problem, etc. Hence, one obtains a set of analytic phases, each one with 
a different number of regular scattering solutions. These phases are separated by 
decay thresholds and, together, form the multichannel coupled OQS. 

Each decay threshold (branching point) can be associated with a non-analytic 
point of the scattering matrix. The coupling of different SM eigenfunctions to the 
same decay channel induces a mixing among them which reflects the nature of the 
branching point. In this way, regular scattering solutions are modifled in the prox- 
imity of each threshold. Such configuration mixing involving all SM eigenfunction of 
the same quantum numbers (angular momentum and parity) can radically change 
the structure of near-threshold SM states. This is the reason for coexistence of 'good 
SM' and 'good CM' states in the spectra of light nuclei. That is also why SM and 
CM schemes can be employed to model wave functions in analytic phases. 

The fiux conservation (unitarity) in a multichannel system implies that the mix- 
ing of SM eigenfunctions changes whenever a new channel opens up. In this sense, the 
eigenfunctions are not immutable characteristics of the system but they strongly de- 
pend on the energy window imposed theoretically and vary with the increasing total 
energy of the system. The role of unitarity and channel coupling for near-threshold 
eigenfunctions has been demonstrated within an OQS formalism for spectroscopic 
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Fig. 2. The schematic picture of couplings in the multichannel representation of ^^C. With the 
increasing excitation energy, subsequent decay channels open up at threshold energies Ei, E2, 
E3, . . . , leading to a complex multichannel network of couplings. When a new channel opens 
up at the threshold Ei, the unitarity imposes the appearance of new channel couplings; hence, 
a modification of all eigenfunctions with energies E < Ei. 



factor^ and radial overlap integrals.'^ 

Figure [2] illustrates the multichannel network of couplings on the example of ^^C. 
For the excitation energies E below the first decay threshold £12^^^^, the regular 
phase contains only elastic scattering solutions. For E > Ei2q^^^ and below the 
second decay threshold at Ei2Q^iiBj^p, the regular phase contains two solutions cou- 
pled by the unitarity condition, and the story goes on at higher excitation energies. 
The structure of any eigenfunction in this multichannel OQS depends on the total 
energy of the system, which in turn determines the environment of decay channels 
accessible for couplings. With increasing excitation energy, the multichannel OQS 
network of couplings becomes increasingly more complex. 

§3. Open quantum system formulation of the nuclear shell model 

In the standard CQS formulation of the nuclear SM, nucleons occupy bound 
single particle orbits of an infinite (harmonic oscillator) potential; hence, they are 
isolated from the environment of scattering states Since the scattering contin- 
uum is not considered, the presence of decay thresholds and exceptional points is 
neglected. Hence, the essential aspect of the nuclear clustering mechanism is to- 
tally absent. In CM, couplings to decay channels are put by hand by introducing a 
posteriori suitable effective cluster degrees of freedom. 

To formulate the nuclear SM for OQSs, two schemes have been proposed. The 
first one, the Gamow Shell Model (GSM),'^^ is the complex-energy CSM based on 
the one-body Berggren ensemble.'^ The GSM, which is conveniently formulated in 
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the Rigged Hilbert Space, ' offers a fully symmetric treatment of bound, resonance, 
and scattering single-particle states. The second scheme, the real-energy 
uses the projection formalism in the Hilbert space. A recent realization of this 
approach is the Shell Model Embedded in the Continuum (SMEC),'^ which offers 
a unified description of the structure and reactions with up to two nucleons in the 
scattering continuum using realistic SM Hamiltonians. In the following, we shall 
employ SMEC to illustrate the mixing of SM eigenstates induced by the coupling to 
the decay channel. 

The detailed description of SMEC can be found in recent reviews.'^''^ The 
Hilbert space is divided into orthogonal subspaces (// = 0, 1, . . . ) corresponding 
to the different number of particles in the scattering continuum. An OQS description 
of 'internal' dynamics, i.e., in Qq, includes couplings to the environment of decay 
channels and is given by the energy-dependent effective Hamiltonian. In the approx- 
imation of one-nucleon scattering continuum, this effective Hamiltonian reads: 

n{E) = Ho + Hi{E) = Ho + V^h{E), (3-1) 

where Hq is the CQS Hamiltonian (the SM Hamiltonian), Vq is the continuum- 
coupling constant, E is the scattering energy, and h{E) is the term describing the 
coupling between localized states (Qo) and the environment of one-nucleon decay 
channels (Qi). The 'external' mixing of two SM eigenstates i and j due to Hi{E) 
consists of the Hermitian principal value integral describing virtual continuum exci- 
tations and the anti-Hermitian residuum that represents the irreversible decay out 
of the internal space Qo- 

An explicit energy dependence of the effective Hamiltonian is the origin of its 
strong non-linearity. Moreover, the continuum-coupling term generates effective 
many-body interactions in the internal space, even if it has two-body character in the 
full space. Therefore, one should not separate a problem of many-body interactions 
in model spaces from a problem which of the two formulations of the many-body 
frameworks is being used: the CQS framework or the OQS framework. 

The biorthogonal eigenstates and of Ti are linear combinations of SM 
eigenstates \ipi): 

where [bji] is an orthogonal matrix. The left \<Pi) and right |^) eigenstates are 
related by the operation of complex conjugation. Each eigenstate of V-iE) is coupled 
to states in neighboring nuclei via a network of reaction channels, either closed or 
open. Contribution of different reaction channels to the total continuum coupling is 
non- uniform and spans over a considerable range of excitation energies 

In general, the continuum coupling lowers binding energies of near-threshold 
CQS eigenstates. The continuum-coupling correlation energy to the SM eigenstate 

Ell;iiE) = mn - H,\^,) ^ Vii^MEm), (3-3) 

depends on the structure of the SM eigenstate and on the nature of decay channels 
involved. 



6 



J. Okolowicz, M. Pioszajczak, and W. Nazarewicz 




Fig. 3. A schematic representation of the continuum-couphng correlation energy correction to the 
CQS eigenstates (SM eigenstates) as a function of the distance from the i = neutron decay 
threshold. The symbols ii and 42 denote two eigenstates within the energy interval characteristic 
of strong coupling to the common decay threshold. 

The correlation energy is peaked at the threshold only for the coupling to a 
neutron decay channel with £ = 0.^^' For higher values and/or for charged parti- 
cles such as protons, deuterons, a-particles, etc., the centrifugal barrier and/or the 
Coulomb barrier shift the maximum of correlation energy above the threshold. 

The continuum coupling correlation energy has a characteristic form of an ap- 
proximately symmetric peak centered either at the threshold for £ = neutron decay 
channel, or above it in the continuum in other cases. The full width at half maximum 
a varies from case to case but a typical value is in the range: 2 MeV< a <3 MeV. 

Figure [3] shows a typical shape of the continuum-coupling correlation energy for 
SM states that are coupled to a neutron decay channel with £ = 0. The maximal 
value of -£'((,oi.r) j ™8'y vary from a few tenths of keV to several MeV, depending on 
the configuration of SM states involved and the nature of the decay channel. The 
relatively small value of the correlation energy is a consequence of the fact that only 
a few SM states of the same quantum numbers enter in the narrow window of strong 
correlations around the decay threshold. 

Even though the continuum-coupling correlation energy is small as compared to 
the binding energy, one should remember that it is of the same order as the pairing 
correlation energy, which profoundly modifies the independent particle motion in 
the atomic nucleus. It is therefore pertinent to ask whether correlations induced by 
the continuum coupling may lead to instabilities of certain SM eigenstates near the 
channel threshold, and whether such couplings represent a collective phenomenon 
involving many SM eigenstates having the same quantum numbers. This question 
will be further addressed in Sec. HI 

For E < (bound system), the eigenvalues Ei{E) of 'H{E) are real. In the 
continuum region, Ei{E) correspond to the poles of the scattering matrix and H 
becomes complex-symmetric. The competition between Hermitian (Hq) and non- 
Hermitian (Hi = VQh{E)) parts of the effective Hamiltonian (|3Tp may lead to 
the coalescence of two eigenvalues, i.e., to the formation of the exceptional point 
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det['H{E;Vo)-£I] = 0, with z^ = 0,1. 



(3-4) 



Single-root solutions of Eq. (|3-4|) correspond to EPs associated with decaying (cap- 
turing) states. For the Hamiltonian (|3-ip . the maximum number of such roots is: 



where n is the number of states of given angular momentum J and parity tt. The 
factor 2 in Eq. (j3-5p comes from the quadratic dependence of 7i on Vq (see Eq. ()3-ip ). 
which yields identical solutions when Vq —Vq. 

The wave function mixing in CQS is closely related to features of exceptional 
points of the complex-extended CQS Hamiltonian. Similarly, to understand the con- 
figuration mixing in an OQS one must analyze the spectrum of EPs of the complex- 
extended OQS Hamiltonian. Since the OQS Hamiltonian (|3Tp is energy dependent, 
the essential information about configuration mixing is contained in the trajectories 
of coalescing eigenvalues Ei^ (E) = Ei^ {E) of the effective Hamiltonian - the so-called 
exceptional threads - for a complex value of the continuum coupling Vq.'^ In the 
following section, we shall describe the mixing of SM eigenvalues close to the proton 
decay threshold in terms of the exceptional threads for the complex-extended SMEC 
Hamiltonian. 

§4. Configuration mixing near the charge particle emission threshold 

Let us discuss the near-threshold configuration mixing and the role of exceptional 
threads on the example of J'^ = 0^, (i = 1, . . . , 4) SM eigenstates of the two-proton 
emitter ^^N^^ coupled to the ^ = proton decay channel leading to the 1/2]^ state 
of ^^F. Our objective is to understand generic aspects of the mixing between different 
realistic CQS eigenstates due to the coupling to the common ^^F(l/2^)+p(£ = 0) 
decay channel rather than reproduce experimental data. 

The SMEC calculations have been carried out in the 0^5/2, lsi/2 SM model 

space. In this space, there are four J'^ = 0^ states in ^^Ne. For Hq we take the 
ZBM Hamiltonian,^^' which correctly describes the configuration mixing around the 
N = Z = 8 shell closure. The residual coupling between Qq and the embedding 
one-proton continuum Qi is generated by the contact force V12 = Vo6{ri — r2). We 
do not consider the two-proton continuum space Q2 so our model ^^Ne is closed to 
two-proton decay. 

The size of the continuum-coupling correction to Hq depends on two parame- 
ters: the continuum coupling strength Vq and the system excitation energy E with 
respect to the threshold {E = 0). The range of relevant Vq values can be deter- 
mined, for example, by fitting decay widths of the lowest states in ^^F. For a ZBM 
Hamiltonian, the experimental decay widths of the ground state l/2f and the first 
excited state 5/2f in ^^F are reproduced by taking Vq = —3500 ± 450 MeV-fm^ and 
Vb = -1100 ± 50 MeV-fm3, respectively. The error bars in Vq reflect experimental 
width uncertainties. 
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4.1. Appearance of the aligned state 

Figure S] shows the predicted energies and widths of the lowest O"*" states of 
-^^Ne. Experimentally, the proton separation energy is very small, Sp ~ 100 keV,l^ 
so all these states are close to the one proton threshold (placed at = in Fig. 
U]). Below the proton threshold, the continuum-coupling correction is Hermitian 
and the eigenvalues of H are real. The lowest energy eigenvalue shows pronounced 
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Fig. 4. Energies Ei (left) and widths Ft (right) of the four lowest J'" = 0^ eigenvalues of the effective 
Hamiltonian (|3-ip for ^®Ne as a function of the proton energy E. The proton threshold at iJ = 
is indicated by a vertical dotted line. The continuum coupling strength is Vo = —1000 MeV fm'^. 

dependence on the continuum coupling in the vicinity of the threshold, as seen in 
Fig. EJ The minimum of Ei is shifted above the threshold by about 0.55 MeV 
due to the Coulomb interaction. The width Fi of all eigenvalues becomes different 
from zero only above the threshold. The maximum of ri(E) is found close to the 
energy corresponding to the minimum of the real part Ei{E). At higher energies, 
ri{E) gradually decreases.'^ As seen in Figs. H] and [5l continuum coupling has 
the strongest effect on the lowest-energy eigenvalue with the largest width and the 
correlation energy rapidly decreases for excited states. 
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Fig. 5. Similar as in Fig. 2] except for the continuum-coupling correlation energy (|3-3[) . 

Figure [6] shows the real part of the off-diagonal matrix elements of the orthog- 
onal transformation [hji] (Eq. ()3-2p ) corresponding to the coupling of 0^, (j = 
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Fig. 6. The real part of the coefficients bji of the orthogonal transformation H3-2|l for Oj" eigenstates 
in ^''Ne. 



1, . . . , 4) SM eigenstates in ^^Ne to the one-proton £ = decay channel. The maxi- 
mum/minimum of hji[E) indicates the region of the strongest mixing at the maxi- 
mum of the continuum-couphng correlation energy (Emax — 0.55 MeV). (The imag- 
inary part of bji{E) exhibits similar behavior.) The large matrix elements involve 
the of state of Ti, which couples strongest to the decay threshold. This single state, 
which we shall call the 'aligned state' in the following, couples strongest to the decay 
channel; hence, it carries many of its characteristics. Aligned states are crucial for 
understanding the mixing of SM eigenstates around the reaction threshold. 

4.2. Exceptional threads 

A complete picture of the near-threshold configuration mixing in CSM many- 
body wave functions can be obtained by investigating properties of exceptional 
threads of the complex-extended CSM Hamiltonian. Figure [7] shows the exceptional 
threads in the complex- Vq plane for the 0^ states of ^^Ne. According to (|3-5p there 
are 12 doubly-degenerate exceptional points for 4 eigenstates. Six of these corre- 
spond to pairs of eigenvalues with the outgoing asymptotics (decaying resonances) 
and six to pairs of eigenvalues with the incoming asymptotics (capturing resonances). 
For E < 0, the effective Hamiltonian (j3Tp is Hermitian; hence, exceptional threads 
for capturing and decaying resonances have identical 9Vo and their real parts have 
opposite signs. In general, threads for 3%Vb • QVq > (KVq ■ ^^^o < 0) correspond to 
decaying (capturing) states. This symmetry does not hold for E > 0. 

Figure EJa) shows six threads in the region of small values of the coupling con- 
stant. Each of those double poles is a result of the coalescence of the lowest-energy 
eigenvalue Of of the effective Hamiltonian with one of the remaining three eigen- 
values. All eigenstates are mutually mixed through their couplings with the aligned 
state of . All exceptional threads relevant to the mixing involve this aligned eigen- 
state. As seen in Fig. [SJ the aligned eigenstate exhausts ~ 80% of the continuum- 
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Fig. 7. Exceptional threads for the 0^ eigenvalues of the SMEC Hamiltonian for ^''Ne in the 
complex-Vb plane. Each point at the exceptional thread corresponds to an exceptional point at 
a definite energy E. The threshold position {E = 0) is indicated for each thread with a circle. 
Because of the invariance of H to Vo — i- — Vb only half of all possible exceptional threads are 
shown in the region of small (a) and very large (b) values of Vq. 

Exceptional threads in (a) formed by the coalescence are marked by solid and short- 

dashed lines for decaying and capturing double-poles, respectively, while the (0^, Oj) and {Of, 
0^) pairs are indicated by long dashed and double-short-dashed, and dotted and dashed-dotted 
lines, respectively. Exceptional threads in (b) do not involve the aligned state 0^. 



coupling correlation energy. This is a generic feature of aligned states, which we also 
found in other systems close to the 1 = proton decay threshold. 

One should notice that the point that is closest to the physical limit of the SMEC 
Hamiltonian (QVq = 0), the turning point, has moved away from the threshold into 
the scattering continuum. This turning point of all exceptional threads involving 
the aligned state appears consistently at the same energy Emax that corresponds to 
a minimum of Ei (see Fig. U]), a maximum of Ecorr in 0^ (see Fig. [5]), and extrema 
oibik{E) (see Fig. E]). 

The remaining six double-poles of the scattering matrix are found at exceedingly 
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large values of complex Vq, outside the range of physically significant values, as seen 
in Fig. El^b). In this case, all exceptional threads for capturing resonances remain in 
the 5?Vo > half-plane also for E > 0. With increasing energy of the system, one 
pair of double poles escapes towards KVq = ±oo, the second pair escapes towards 
QVo = iboo, and the third one remains confined in a relatively narrow range of Vq. 

4.3. Maximum continuum coupling point 

As discussed above, the continuum-mixing among 0^^ SM eigenstates of ^^Ne is 
strongest at the turning point at Emax — 0.55 MeV above the £ = proton-decay 
threshold. In general, the scattering energy corresponding to the turning point 
depends on the height of the effective (Coulomb+centrifugal) barrier, structure of 
SM wave functions, and the channel itself. In the studied case of ^ = decay, the 
centrifugal barrier is absent. The height of the Coulomb barrier for Zt = Z — 1 = 9 
in ^^Ne is Vc* = 1.69 ± 0.01 MeV. This estimate has been obtained by considering 
a uniformly charged sphere of radius Rq = 1.27{A — 1)^/^ fm. The error bar on Vc 
corresponds to the variation C/q = 52 it 5 MeV of the depth of the Woods-Saxon 
(WS) potential of radius Rq and diffuseness d = 0.67 fm. 

To see the dependence of Emax on Vc, we investigate the mixing of two 0+ SM 
states in ^^S due to the coupling to the one-proton decay channel ^'^P(l/2^)+p(^ = 
0). In this case, there are two exceptional threads associated with decaying and 
capturing resonances. SMEC calculations are performed in the lsi/2> 0^5/2, 0^3/2 
SM model space using the USDB interaction.!^ 
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Fig. 8. The real part of &12 in {Zt = 15, solid line). The dashed line shows results with the 
Coulomb potential reduced by a factor of three {Zt = 5). The turning points are marked by 
arrows. 

Figure [8] shows the energy dependence of the real part of the mixing matrix 
element buiE). For ^^s, the height of the Coulomb barrier is Vc = 2.76 ±0.03 MeV. 
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To assess the impact of the Coulomb barrier, we also carried out calculations with 
the reduced Coulomb potential having Vc = 0.82 it 0.007 {Zt = 5). One can see that 
Emax strongly depends on Vc- The turning point of exceptional threads in ^^S is 
1.25 MeV and 0.17 MeV for standard and reduced Coulomb potentials, respectively. 
Notice that the magnitude of mixing, as given by the amplitude 1 612 (E) \ , is strongly 
reduced with increasing charge Zt. This means that the imprint of the branching 
point at a charge particle decay threshold on regular solutions of a scattering problem 
is shielded by the Coulomb barrier; hence, it is reduced in heavier nuclei. This generic 
feature, with far-reaching consequences for the manifestation of nuclear clustering, 
does not depend on the nature of the decay channel. 




Fig. 9. Exceptional threads for the 0^ states in ^^S in the complex- Vo plane for the standard (solid 
line) and reduced (dashed line) Coulomb barrier. 



Finally, Fig. [9] shows the exceptional threads in the complex- Vq plane for the 
0"^ states in -^^S. We show here only those exceptional threads which are associated 
with decaying resonances. It is seen that the turning point moves further away from 
the physical limit at 9Vo=0 with increasing charge; this nicely illustrates the find- 
ings of Fig. [HI i.e., the continuum-coupling induced configuration mixing diminishes 
with Zf. This generic effect due to the Coulomb barrier not only suppresses proton 
halos in heavier nuclei, but also prevents a clustering near a emission thresholds in 
heavy systems. Conventionally, the latter effect has been associated with the disrup- 
tive influence of the nuclear spin-orbit interaction on a-clustering (see Ref.'^ and 
references cited therein), i.e., with a specific feature of the nuclear potential. 
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§5. Conclusions and perspectives 

The phenomenon of nuclear chistcring has been a key problem in low-energy 
nuclear structure studies. In this work, we have applied the OQS framework of 
SMEC to study the behavior of SM eigenstates in the vicinity of the i = proton- 
decay threshold. Proximity of the branching point induces the collective mixing of 
SM eigenstates, in which the essential role is played by a single aligned eigenstate of 
the OQS Hamiltonian. The aligned state, in which most of the continuum-coupling 
correlation energy is concentrated, is present in all exceptional threads close to the 
physical limit of the OQS Hamiltonian (S>Vb=0) for realistic values of the real part of 
the coupling constant 3?Vb. Above the decay threshold, this state becomes a broad 
resonance. 

The collective mixing of SM eigenstates via the aligned state - an imprint of the 
decay channel - is the essence of the clustering phenomenon. In the studied exam- 
ples of the coupling to i = proton decay channel, an aligned state can cause an 
instability of a near-threshold SM eigenstate if its energy is close to the turning point 
of associated exceptional threads. We thus conjecture that clustering present in the 
vicinity of any charged-particle emission threshold is a consequence of a mixing medi- 
ated by an aligned state involving all SM eigenstates with the same quantum numbers 
J'^. Hence, the presence of charged-cluster states near the respective charged-cluster 
emission threshold is a signature of a profound change of the near-threshold SM wave 
function and the direct manifestation of continuum-coupling correlations. 

In order to sec an imprint of the aligned state in a near-threshold SM wave func- 
tion, the distance of the turning point from a real axis 3'Vo=0 should be small. This 
distance, controlled mainly l)y the Coulomb barrier penetrability, increases rapidly 
with increasing proton number. This feature, which is generic for all configura- 
tion mixing situations involving the i = charged-particle (e.g. proton, deuteron, 
a-particle) emission, reduces clustering correlations in heavier nuclei. 

Favorable conditions for the appearance of charged-cluster configurations are 
above the charged-cluster emission threshold. On the other hand, neutral-cluster 
configurations are expected to show up primarily below the threshold - due to the 
rapid growth of the decay width with energy. Here, spectacular examples are one- 
and two-neutron halos in light nuclei. 

In summary, we have shown that the many-body correlations due to continuum 
coupling can profoundly impact the nature of states near reaction thresholds. While 
much work is still needed to fully understand the phenomenon of nuclear clustering, 
the recent progress in microscopic nuclear theory of OQSs, and insights provided by 
CSM, make us optimistic about the solution to this fascinating puzzle. 
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